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We review recent developments in charm physics, focusing on the physics of charmed 
mesons. We discuss charm spectroscopy, decay constants, as well as searches for new 
physics with charmed mesons. We discuss — mixing and CP-violation in charm 
decays. We also present the modified Nelson plot of charm mixing predictions. 



1 Introduction 

Charm physics plays a unique dual role in the modern investigations of flavor physics. Charm 
decay and production experiments provide valuable checks and supporting measurements for 
studies of CP-violation in measurements of CKM parameters in b-physics, as well as outstanding 
opportunities for searches for new physics. Historically, many methods of heavy quark physics 
have been first tested in charmed hadrons. The fact that a b-quark mainly decays into a charm 
quark makes charm physics an integral part of any b-physics program. In many cases, direct 
measurements of charm decay parameters directly affect the studies of fundamental electroweak 
physics in B decays p. 

This year brought several interesting developments in some seemingly well-understood sec- 
tors of charm physics, such as meson spectroscopy. Here I shall discuss theoretical implications 
of these and other results. The experimental status of charmed meson spectroscopy was dis- 
cussed in R. Chistov's talk . Resent results in the measurements of charmed meson formfac- 
tors, lifetimes and — mixing parameters were discussed by W. Johns |2| and G. Boca 



2 Spectroscopy 

Meson spectroscopy is an important laboratory for understanding quark confinement. Mesons 
containing one heavy quark can provide valuable information about the structure of the QCD 
Lagrangian, as spectroscopic considerations simplify significantly in the limit of infinitely heavy 
quark, rriQ/A —>■ oo, where A represents a typical scale of hadronic interactions. While charm 
quark hardly satisfies this conditions, it is nethertheless useful to apply these considerations 
to the charmed quark systems. In this limit the heavy quark spin Sq decouples, so the total 
angular momentum of the light degrees of freedom Jf becomes a "good" quantum number. 
Since parity of a meson can be obtained by knowing the angular momentum quantum number 
/ as (— l)'"*"^, this leads to an important prediction of heavy quark symmetry: the appearance 
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of heavy meson states in the form of degenerate parity doublets classified by the total angular 
momentum of the light degrees of freedom (see Table (HJ), 

s' = Jf±l. (1) 

This mass degeneracy is lifted with the inclusion of subleading I/tuq corrections. This useful 
picture is built into many quark-model descriptions of heavy meson spectra. The resulting 
models have been very successful in explaining the spectrum of negative-parity scalar and 
vector Jf = l/2~ and positive-parity vector and tensor Jf = 3/2"'" states. A narrow resonance 
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1/2 


1/2 
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3/2 


5/2 


S 


0,1 


0,1 


1,2 


1,2 


2,3 



Table 1: Total angular momentum assignments for heavy-light mesons 

in -D^vr" was recently reported by BaBar and confirmed by the CLEO ^ and Belle |Zj 
collaborations. Its decay patterns suggest a quark-model 0"*" classification, which would identify 
it with the positive-parity Jf = 1/2+ p-wave state. As in the D meson system, p-wave states for 
the Df system are expected, and two narrow states, Dsi(2536) and Ds2(2573) were discovered 
by ARGUS and CLEO, respectively JHj. In analogy to the D system, two broad states are also 
expected. 

The mass of the new state 2317.6 ± 1.3 MeV appears surprisingly low and its width appears 
to be too small for quark model practitioners. In fact, this state appears below DK threshold, 
closing off the most natural decay channel for this state. This forces it to decay mainly via an 
isospin-violating transition into the Dfir^ final state which makes its width quite narrow. Its 
mass disagrees with most predictions of quark models UHl IHl IIS CSl El available prior to its 
observation. For example, a mass of 2487 MeV is obtained in the potential model calculation 
by Eichten and Di Pierro [lOj. Quenched lattice calculations also seem to favor larger values 
of the mass of this state ^H] (see, however, jEHl)- This led to a lively discussion of the possible 



Reference 


0+ mass 


1"'" mass 


Ebert et al (98) 14 
Godfrey-Isgur (85) P 
DiPierro-Eichten (01) ^ 
Gupta- Johnson (95) [TT] 
Zeng et al (95) ^ 
Experiment 


2.51 GeV 

2.48 GeV 

2.49 GeV 
2.38 GeV 
2.38 GeV 
2.317 GeV 


2.57 GeV 
2.55 GeV 
2.54 GeV 
2.52 GeV 
2.51 GeV 
2.463 GeV 



Table 2: Theoretical predictions for masses of 0"*" and 1+ Dg states 

non-gg nature of this state |18| IT^ 1201 j^. A possibility of a state that is an admixture of a 
four-quark state and a qq states was discussed in j2Zl- In addition, a second narrow state is 
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observed in D*+7r° at a mass near 2460 MeV [HI EI- This state would naturally be identified as 
a spin 1 positive parity p-wave meson. However, its mass also appears too low for the potential 
model expectations (e.g. 2605 MeV jTHl EH], see also Table 121). Its radiative decays to the 
ground state Dg meson were observed with 

r(D..(2460) ^ ^ 44^010 .2) 

r(D.,(2460) - Z^.ttO) U.44±0.iO, {I) 

while doubly charged states were not observed in Dfn^ channels. Finally, Dgj states were also 
observed in B-decays B DD^*] . 

The low values of the masses for these states, however, do not signal a breakdown of quark- 
model descriptions of the heavy meson spectrum, as it is difficult to assess the accuracy of these 
predictions, especially in the charm sector. Many authors make use of the non-relativistic nature 
of the charm quark, taking into account l/rric corrections only. For the O"*" state, quark model 
predictions range from the values of 2387 — 2395 MeV fUEj (still above the DK threshold) on 
the low end of the spectrum to 2508 MeV [ll] on the high end. Since the described phenomena 
are highly non-perturbative, one should be careful before making a judgment on the nature of 
a given state based solely on the prediction of a given quark model. For example, as discussed 
above, in the mc — oo limit the 0"*" and 1"*" states are expected to become degenerate in mass, 
mo+ , mi+ — i> M. This can be emulated in quark models by neglecting heavy-quark symmetry- 
violating l/rric corrections. Yet, different quark models predict very different behavior in this 
"heavy-quark limit": for instance, one potential model fTT] predicts that the mass M of the 
(O"*", l"*") multiplet will decrease to approximately 2382 MeV (which is less than the mass of the 
0"*" state predicted in this model with the full potential), while in a QCD string model jTHI it 
is expected to increase up to 2500 MeV (which is much greater than the mass of the O"*" state 
predicted in this model with the full potential). In addition, quark models, modified to include 
chiral symmetry constraints, generally predicted lower values of mass splitting between (0~, 1~) 
and (0^, 1^) multiplets, of the order of 200 — 300 MeV |17j. In addition, one has to remember 
that most of the unusial details about these states, such as narrowness of their decay widths, 
simply follows from the fact that the mass of that state is smaller than the D^*^K threshold. It 
is then only the fact that the new state appears below DK threshold and is almost degenerate 
with a non-strange O"*" p-wave D state [221 curious and deserves an investigation, although 
could be purely accidental. 

A combination of experimental measurements described above can shed some light onto the 
nature of these states. For instance, molecular-type explanation of the low masses of these 
states implies the existence of the doubly-charged states, which were not observed. On the 
contrary, a possible disagreement of the observed branching ratios of B decays into these states 
with calculations of their branching ratios in naive factorization could favor molecular nature 
of these states [21] • But it could as well signal a breakdown of naive factorization in B decays 
into the pair of 0~,0+(l^) open-charm states, which was never really tested, or simply refiect 
our ignorance of the decay constants of positive parity mesons [20]. The observation of the 
radiative decay of Dsj{24Q0) favors qq (or maybe gg-four-quark-state admixture) explanation 
of the nature of these mesons. 
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3 Decay constants and B-physics experiments. 



Since m^, rric ^ ^qcd, both charm and bottom quarks can be regarded as heavy quarks ^. 
Naturally, heavy quark symmetry relates observables in B and D transitions. As an example, 
let us consider measurements in the charm sector affect determinations of the CKM matrix 
elements relevant to top quark in — mixing. 

A mass difference of mass eigenstates in B^ — B^ system can be written as 



Am^ =C \af\^l) 



-6/23 



1 + 
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(5) 
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{b:\o{M), 



(3) 



where C 



GpM^,{Vtb*Vtd) 'r]BmBSo{xt)/ (47r^) (see Ref. ^ for complete definitions of the 



parameters in this expression ). The largest uncertainty of about 30% in the theoretical calcu- 
lation is introduced by the poorly known hadronic matrix element A = (5^|0(/x)|i?°). Eval- 
uation of this matrix element is a genuine non-perturbative task, which can be approached 
with several different techniques. The simplest approach ("factorization") reduces the ma- 
trix element A to the product of matrix elements measured in leptonic B decays A-^ = (8/3) 
{B^\bL'y„dL\0){0\bL'-f'^dL\B'^) = (2/3)/^m^, where we employed the definition of the decay con- 
stant fs, 

{0\hi,dL\B%p))=zpjB/2. (4) 

A deviation from the factorization ansatz is usually described by the parameter Bb^ defined 
as y4 = Bb^A-^; in factorization Bb^ = 1- Similar considerations lead to an introduction of the 
parameter Bb^ defined for mixing of Bg mesons. It is important to note that the parameters 
BBg depend on the chosen renormalization scale and scheme. It is convenient to introduce 
renormalization-group invariant parameters Bb„ 
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We provide averages oi Bb^, as well as the ratio BbJBb^ from the review [32] as well as from 
two more recent evaluations |33tlM] in Table El Thus, at least naively, one can determine CKM 
matrix element Vtd by measuring /b and Am^ and computing Bb^- 

This direct approach, however, meets several difficulties. First, leptonic decay constant Jb 
can in principle be extracted from leptonic decays of charged B mesons. The corresponding 
decay width is 



T{B ^ lu) = -^l 



2 IT/ |2 2 

B \ Vub\ miniB 



1 



mf 



m 



(6) 



B, 



This width is seen to be quite small due to the smallness of the CKM factor \Vub\ and helicity 
suppression factor of rrif. In addition, experimental difficulties are also expected due to the 
backgrounds stemming from the presence of a neutrino in the final state. 



"'^This approximation obviously works better for bottom than for charm quarks. 
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Method (reference) 




BbJBb^ 


Lattice, '03 [32j 
QCDSR, '03 
Lattice, '03 ^ 
QCDSR, '03 jSU 


1.34(12) 
1.67 ±0.23 
1.277(88)(li) 
1.60 ±0.03 


1.00(3) 
^ 1 

1.017(16)(lf^) 
^ 1 



Table 3: Renormalizat ion-group independent B-parameters. 



Second, computation of Bb^ is quite difficult and requires the use of non-perturbative tech- 
niques such as lattice or QCD Sum Rules. Current uncertainties in the determinations of 
and Bb„ are quite large. It turns out that evaluation of the ratio of 



Am^ _ rriBa 
Am, 



BbJ. 



B, 



V Bbs Jbs 



Vtd ' 
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(7) 



is favored by lattice community, as many systematic errors cancel in this ratio. This gives a 
ratio of which provides a non-trivial constraint on CKM parameters in the p — rj plane. 

Instead, one can make use of ample statistics available in charm production experiments, 
as heavy quark and S't/(3) flavor symmetries relate the ratio of charm decay constants fos/fD 
to beauty decay constants fss/ Ib 



Ibs I fh 



D 



1 ± 0{ms) X 0(l/m6 - l/mj. 



Note that SU(3)-violating corrections can also be evaluated in chiral perturbation theory [33]. 
One still needs to rely on the theoretical determination of Bb^- 

Similar techniques of relating B and D decays can also be used to extract other CKM 
matrix elements, hke Vub studies of lifetime patterns of heavy hadrons [HTj, and tuning 
lattice QCD calculations 



4 Charm mixing and CP violation 

One of the important areas of modern phenomenology where charm decays play an important 
role is the indirect search for physics beyond the Standard Model. Indeed, large statistics 
usually available in charm physics experiment makes it possible to probe small effects that 
might be generated by the presence of new physics particles and interactions. A program of 
searches for new physics in charm is complimentary to the corresponding programs in bottom 
or strange systems. This is in part due to the fact loop-dominated processes such as — 
mixing or flavor-changing neutral current (FCNC) decays are sensitive to the dynamics of ultra- 
heavy down-type particles. Also, in many dynamical models, including the Standard Model, 
the effects in s, c, and h systems are correlated. 

The low energy effect of new physics particles can be naturally written in terms of a series 
of local operators of increasing dimension generating AC = 1 (decays) or AC = 2 (mixing) 
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transitions. For — mixing these operators, as well as the one loop Standard Model effects, 
generate contributions to the effective operators that change D° state into D*^ state leading to 
the mass eigenstates 

|Di)=p|D°)±g|:D°), (9) 

where the complex parameters p and q are obtained from diagonalizing the D° — mass 
matrix. The mass and width splittings between these eigenstates are given by 

1712 -mi Fs-Ti 

y^^^- (10) 

It is known experimentally that — D° mixing proceeds extremely slowly, which in the 
Standard Model is usually attributed to the absence of superheavy quarks destroying GIM 
cancellations [39j. 

It is instructive to see how new physics can affect charm mixing. Since the lifetime difference 
y is constructed from the decays of D into physical states, it should be dominated by the 
Standard Model contributions, unless new physics significantly modifies AC = 1 interactions. 
On the contrary, the mass difference x can receive contributions from all energy scales. Thus, 
it is usually conjectured that new physics can significantly modify x leading to the inequality 
X ^ y ^. The same considerations apply to FCNC decays as well, where new physics could 
possibly contribute to the decay rates oi D ^ X^'j, D XJ^l^ (with X„ being exclusive or 
inclusive final state) as well as other observables One technical problem here is that in 
the standard model these decays are overwhelmingly dominated by long-distance effects, which 
makes them extremely difficult to predict model-independently. This problem can be turned 
into a virtue |^ . 

Another possible manifestation of new physics interactions in the charm system is associated 
with the observation of (large) CP- violation. This is due to the fact that all quarks that build up 
the hadronic states in weak decays of charm mesons belong to the first two generations. Since 
2x2 Cabbibo quark mixing matrix is real, no CP- violation is possible in the dominant tree-level 
diagrams that describe the decay amplitudes. In the Standard Model CP-violating amplitudes 
can be introduced by including penguin or box operators induced by virtual 6-quarks. However, 
their contributions are strongly suppressed by the small combination of CKM matrix elements 
VcbV*fj. It is thus widely believed that the observation of (large) CP violation in charm decays or 
mixing would be an unambiguous sign for new physics. This fact makes charm decays a valuable 
tool in searching for new physics, since the statistics available in charm physics experiment is 
usually quite large. 

As in B-physics, CP-violating contributions in charm can be generally classified by three 
different categories: (I) CP violation in the decay amplitudes. This type of CP violation 
occurs when the absolute value of the decay amplitude for D to decay to a final state / (Af) 
is different from the one of corresponding CP-conjugated amplitude ("direct CP- violation" ) ; 
(II) CP violation in — mixing matrix. This type of CP violation is manifest when 
Ki = \p/q.? = (2Mi2 - iri2)/(2M*2 - iVl^) ^ 1; and (III) CP violation in the interference 
of decays with and without mixing. This type of CP violation is possible for a subset of final 
states to which both and can decay. 



^This signal for new physics is lost if a relatively large y, of the order of a percent, is observed |40[ I41| . 
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For a given final state /, CP violating contributions can be summarized in the parameter 



qAf_ 



Rr, 



A, 



where A^ and Af are the amplitudes for — > / and — > / transitions respectively and 8 is 
the strong phase difference between A^ and A^. Here represents the convention-independent 
weak phase difference between the ratio of decay amplitudes and the mixing matrix. 

Presently, experimental information about the — mixing parameters x and y comes 
from the time- dependent analyses that can roughly be divided into two categories. First, more 
traditional studies look at the time dependence of D ^ f decays, where / is the final state 
that can be used to tag the flavor of the decayed meson. The most popular is the non-leptonic 
doubly Cabibbo suppressed decay — > K^n'. Time-dependent studies allow one to separate 
the DCSD from the mixing contribution K^-' 



TT 



-Ft I 



A 



K- 



cos< 



TT 



X smi 



)n + ^{y' + x'){nf 



(12) 



where R is the ratio of DCS and Cabibbo favored (CF) decay rates. Since x and y are small, 
the best constraint comes from the linear terms in t that are also linear in x and y. A direct 
extraction of x and y from Eq. (fT^ is not possible due to unknown relative strong phase 6d of 
DCS and CF amplitudes |l4j, as x' = xcosSd + ysmSc, y' = ycosSn — xsinSn- This phase 
can be measured independently. The corresponding formula can also be written |40J for 
decay with x' —a;' and Rm R^- 

Second, mixing can be measured by comparing the lifetimes extracted from the analysis 
of D decays into the CP-even and CP-odd final states. This study is also sensitive to a linear 
function of y via 



t{D K- 



t{D K+K- 



y cos I 



xsm< 



R^ — 1 



(13) 



Time-integrated studies of the semileptonic transitions are sensitive to the quadratic form x^+y'^ 
and at the moment are not competitive with the analyses discussed above. 

The construction of new tau-charm factories CLEO-c and BES-III will introduce new time- 
independent methods that are sensitive to a linear function of y. One can again use the fact 
that heavy meson pairs produced in the decays of heavy quarkonium resonances have the useful 
property that the two mesons are in the CP-correlated states |45j . 

By tagging one of the mesons as a CP eigenstate, a lifetime difference may be determined 
by measuring the leptonic branching ratio of the other meson. Its semileptonic width should be 
independent of the CP quantum number since it is flavor specific, yet its branching ratio will 
be inversely proportional to the total width of that meson. Since we know whether this D{k2) 
state is tagged as a (CP-eigenstate) D± from the decay of D{ki) to a final state 5*0- of definite 
CP-parity o" = ±, we can easily determine y in terms of the semileptonic branching ratios of 
D±. This can be expressed simply by introducing the ratio 



Rt 



T[^L -^{H^ S,){H ^ X)] Br{HO ^ Xlu) 



(14) 
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where X in H ^ X stands for an inclusive set of all final states. A deviation from = 1 
implies a lifetime difference. Keeping only the leading (linear) contributions due to mixing, y 
can be extracted from this experimentally obtained quantity, 

ycos<P={-l)'a^^. (15) 

The current experimental upper bounds on x and y are on the order of a few times 10^^, 
and are expected to improve significantly in the coming years. To regard a future discovery of 
nonzero x or ?/ as a signal for new physics, we would need high confidence that the Standard 
Model predictions lie well below the present limits. As was recently shown |1T] , in the Standard 
Model, X and y are generated only at second order in SU(3)i? breaking, 

X , 2/ ~ sin^ Oc X [SU{3) breaking]^ , (16) 

where 6c is the Cabibbo angle. Therefore, predicting the Standard Model values of x and y 
depends crucially on estimating the size of SU(3)ir breaking. Although y is expected to be 
determined by the Standard Model processes, its value nevertheless affects significantly the 
sensitivity to new physics of experimental analyses of D mixing W0\. 

Theoretical predictions of x and y span several orders of magnitude. The predictions ob- 
tained in the framework of the Standard Model are not exception, as evidenced from Fig. 
Roughly, there are two approaches, neither of which give very reliable results because rric is in 
some sense intermediate between heavy and light. The "inclusive" approach is based on the 
operator product expansion (OPE). In the rric ^ A limit, where yl is a scale characteristic of 
the strong interactions, AM and AF can be expanded in terms of matrix elements of local 
operators IIHI- Such calculations yield x,y < 10~^. The use of the OPE relies on local quark- 
hadron duality, and on A /rric being small enough to allow a truncation of the series after the 
first few terms. The charm mass may not be large enough for these to be good approximations, 
especially for nonleptonic D decays. An observation of y of order 10~^ could be ascribed to 
a breakdown of the OPE or of duality, but such a large value of y is certainly not a generic 
prediction of OPE analyses. The "exclusive" approach sums over intermediate hadronic states, 
which may be modeled or fit to experimental data jlHI- Since there are cancellations between 
states within a given SU{3) multiplet, one needs to know the contribution of each state with 
high precision. However, the D is not light enough that its decays are dominated by a few final 
states. In the absence of sufficiently precise data on many decay rates and on strong phases, 
one is forced to use some assumptions. While most studies find x,y < 10~^, Refs. ^\ obtain 
X and y at the 10~^ level by arguing that SU(3)i? violation is of order unity, but the source 
of the large SU(3)i? breaking is not made explicit. It was also shown that phase space effects 
alone provide enough SU(3)ir violation to induce y ~ 10~^ jUj. Large effects in y appear 
for decays close to D threshold, where an analytic expansion in SU(3)i? violation is no longer 
possible. Thus, theoretical calculations of x and y are quite uncertain, and the values near the 

■^Compilation of the — Z?" mixng predictions is known as the Nelson plot ^Sj. In order to obtain a 
compilation of the Standard Model (Fig.QJ and new physics (Fig.[21) predictions for charm mixing, we updated 
and corrected Ref. 0^1 to remove double counting of predictions. We also separated the Standard Model and 
new physics predictions into two separate plots |47| . 
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current experimental bounds cannot be ruled out. Therefore, it will be difficult to find a clear 
indication of physics beyond the Standard Model in — mixing measurements alone. The 
only robust potential signal of new physics in charm system at this stage is CP violation. 
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Figure 1: Standard Model predictions for (open triangles) and \y\ (open squares). Horizontal 
line references are tabulated in Table |31 
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Table 4: Theoretical predictions for mixing parameters (Standard Model). The notation "±" 
indicates the range of predictions. 
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